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ABSTRACT 

Extracting sources with low signal-to-noise from maps with structured background is 
a non-trivial task which has become important in studying the faint end of the submil- 
limetre number counts. In this article we study source extraction from submillimetre 
jiggle-maps from the Submillimetre Common-User Bolometer Array (SCUBA) using 
the Mexican Hat Wavelet (MHW), an isotropic wavelet technique. As a case study we 
use a large (ll.Sarcmin^) jiggle-map of the galaxy cluster Abell2218, with a 850-/im 
la r.m.s. sensitivity of 0.6-1 mJy. We show via simulations that MHW is a powerful 
tool for reliable extraction of low signal-to- noise sources from SCUBA jiggle-maps and 
nine sources are detected in the A2218 850- /j,m image. Three of these sources are identi- 
fied as images of a single background source with an unlensed flux of 0.8 mJy. Further, 
two single-imaged sources also have unlensed fluxes < 2 mJy, below the blank-field 
confusion limit. In this ultradeep map, the individual sources detected resolve nearly 
all of the extragalactic background light at 850 //m, and the deep data allow to put an 
upper limit of 44 sources per arcmin^ to 0.2 mJy at 850 //m. 
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1 INTRODUCTION 

The arrival of th e Submillimetr e Common-User Bolometer 
Array fSCUBA: [Holland et al.l TlQQQ') on the James Clerk 
Maxwell Telescope (JCMT) heralded a new era in galaxy 
evolution studies. The p opulations of submi ll imetre-bright , 
high redshift galaxies llSmail et al.l Il997al : iHughes et alJ 
ll998l:lBarger et alJll99^ have been studied using three dif- 
ferent observational strategies; gravitational lens surveys, 
where pointings are made in the regions of well known clus- 
ter lenses hence benefiting from the cluster magnification, 
blank field surveys, generally selected to overlap with areas 
already observed at other wavelengths, and pointed photom- 
etry of targets selected in other bands. 

The blank fi eld surveys iBarger et al.l 
Hughes et aP Il998l: lEales et all I2OO0I: IScott et al.l 



199c ; 



2002 : 



Borvs et alJl2003HWebb et al.ll2003^ have produced counts 
at fluxes brighter than Ssso ~ 2 mJy. This flux limit is 



caused by the source confusion, which for the 15" angular 
resolution of the JCMT at 850 /im appears to become a 
problem at ab out 2 mJy, effectively limiting the depths of 
these surveys JCondonlll974HBlain et alJll998l:lBarger et alJ 
I1999I: lEales et al.l I2OO0I: iHogd I2OOII) . To this flux limit 
35-45% fe.g.■ lBon^^^^lJl2003^ of the submil l imetre extra- 
galac tic background light llPuget et ai]ll996l : iFixsen et alJ 
Il998h is resolved. 

In comparison, the cluster gravitational lens surveys 
(' Sm afl et alJl997al . l20oi : ICowie et alJ2002l: IChaoman et alJ 

■2002i : Knudsen et al., in prep.), probe sources fainter than 
the blank field confusion limit. As the lensing transformation 
is a mapping from source plane to image plane, it changes 
the number density of observed sources; Both the flux of a 
given source and the area of the region surveyed are mag- 
nified. As a result the confusion limit is moved to a fainter 
fiux level, thereby enabling us to reliably study the fainter 
sub-mm galaxy population, which would otherwise not be 



2 K. K. Knudsen, V. E. Barnard et al. 



accessible with the current telescopes and instrumentation. 
The magnification can be very large near critical lines, where 
we e xpect to see multiple images of the background sources 
llKneib et al. 1996). The typical size of the critical lines is 
20-30" in radius for massive clusters and hence SCUBA's 
2.3' diameter field-of-view is very well-matched to the high 
magnification region of typical cluster lenses. 

This faint population, which can currently only be 
probed via lensing, is crucial to study as it is an important 
contributor to the extraealac t ic submillimetre b a ckground 
light (e.g., iBlainI Il997l 119991 iBorvs et alj l2003t iKnudsenl 
|2004) . While many tens of submillimetre galaxies with 
fiuxes >4mJy have already been studied in great detail 
probing their redshift distribution, ga s content etc. (e.g., 
IChapman et all2005l : lGreve et al Eooi), very little is known 
about the faint population as it yet has been possible only 
to identify and st u dy a ha ndf ul of o bjects in detail (e.g., 
iKneib et al ] |2004al: iBorvs et al.ll2004^ . 

In this paper we present one of the deepest ever SCUBA 
maps of a massive galaxy clu ster: Abell 2218, a rich galaxy 
cluster at redshift z = 0.175 jAbefi et al.lll989^ and acts as 
a powerful gravitational lens both distortin g and magnify- 
ing the distant galaxies iKneib et alJll99^ . It has already 
proven to be a powerful tool in studying high redshift, faint 
galaxies. Exploiting the gravitational lensing, several back- 
ground galaxies have been found, including the lensed z — 
5.57 galaxy of Ell is et al, ( 20QL) . Probably the most unusual 
galaxies have been a multiply-i maged submillimetr e galaxy 
seen as three distinct images teneib et alJl2004a^ ■ which 
have ah eady been studied in great detail at other wave- 
lengt hs (iKneib et alJl2005l: ISheth ct al."2004': 'Garrett et al.! 
| 2005^. and a red shift z ~ 7 galaxy IKneib et al.. ■2004bi: 
Eeami et aDl2005t) . 

To analyse this new SCUBA map, we used a novel tech- 
nique based on the Mexican Hat Wavelet (MHW, so-called 
because the mother wavelet is like a so mbrero in 3-d space ) 
technique and rout ines de veloped bv [Cavon et all i200d) . 
IVielva et al] ll2001allbl) and lVielval (120031) for use with antic- 
ipated Planck Surveyor^ data. In the analysis of the Planck 
data MHW will be used for removing contaminating fore- 
ground extragalactic point sources and leave the Cosmic Mi- 
crowave Background (CMB) signal. However, in this paper 
we are interested in the point sources themselves as the end 
product. The key advantage of the MHW technique is that 
without needing to characterise the background and noise in 
an image, Gaussian shaped beams can be picked out from 
a variety of backgrounds with accurate positions and fluxes. 
This makes the technique ideal for submillimetre images 
where sources are significant in size with respect to the field- 
of-view, and the noise is often relativel y high, variable across 
the field and hard to model accuratelv. lBarnard et all (l2004l) 
successfully applied MHW as a source extraction tool on 
SCUBA scan-maps of Galactic regions and placed valuable 
upper limits for the source counts of submillimetre galaxies 
at 50 and 100 mjy. 

This paper is Paper I in a series of two papers and is an 
introduction to the use of the MHW technique in SCUBA 
jiggle-maps followed by its application on the deep Abell 
2218 SCUBA map. The SCUBA data presented here include 
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m ore observations tha n were published for the same cluster 
in lKneib et al.l (j20043). Paper II (Knudsen et al., in prepara- 
tion) presents the multi-wavelength identification of the un- 
derlying galaxies giving rise to the submillimetre emission. 
The present paper is organized as follows: Section 2 covers 
the observations and the data reduction. In Section 3 we dis- 
cuss in details the MHW technique. The detailed analysis of 
the Abell 2218 sub-mm maps, the detected sources and de- 
rived number counts are discussed in section 4. Conclusions 
and prospects are summarised in section 5. In a further pa- 
per (Knudsen et al., in prep.) we will describe the results of 
applying the MHW technique to a wider survey of cluster 
lenses (the Leiden-SCUBA Lens Survey). 

Throughout we will assume an f2 = 0.3, A — 0.7 cos- 
mology with Ho = 70kms~^ Mpc~^. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

We obtained observations of the cluster of galaxies A2218 
at 850 ^im and 450 /im with SCUBA during March 1998, 
August 2000, January 2001 and January 2002. SCUBA is 
a submillimetre mapping instrument operating at 850 /im 
and 450 ^im simuhaneously (Holland et al. 1999). SCUBA 
consists of two arrays of bolometers having 37 elements at 
850 nm and 91 elements at 450 /im. The field-of-view on 
the sky, which is approximately the same for both arrays, is 
roughly circular with a diameter of 2.3'. The observations 
were carried out in jiggle mode, where the secondary mir- 
ror follows a 64-point jiggle pattern in order to fully sample 
the beam at both operating wavelengths. In order to cover a 
larger sky area four pointings were made, though with a big 
overlap region. As the region with large gravitational lens- 
ing magnification is extended due to the double-peaked mass 
distribution of A2218, four pointings were obtained. The off- 
sets between the central position of each pointing were about 
0.5 arcmin. The primary sky subtraction is achieved by chop- 
ping the secondary mirror. The chop configuration used for 
most of the time was a chop throw of 35" with a position an- 
gle fixed in right ascension (RA), though larger chop throws 
and other position angles were also used for a limited part 
of the observations. The resulting beam, weighted by time 
spent using the different chop throws, is shown in Figure 
As a result of chopping to either side the beam pattern has 
a central positive peak with a negative dip on each side in 
the chop direction with amplitudes of half the peak fiux. 

During the observations the pointing was checked every 
hour by observing nearby bright blazars with known posi- 
tions. The noise level of the bolometer arrays was checked 
at least twice during an observing shift, and the atmo- 
spheric opacities at 850 /im and 450 /im, rgso^mand r450pm, 
were determined by a JCMT skydip observation every two- 
three hours and supplemented with the r225GHz data from 
the neighboring Caltech Submillimeter Observatory (CSO). 
The zenith opacity was 0.15 < rgso^m < 0.40 (correspond- 
ing to 0.037 < r225GHz <0.1), although data taken when 
TgsoMm > 0.32 were not included in the 450-/im map since it 
is much more sensitive to poor weather. Calibrators with 
accurately known flux densities were observed every two 
to three hours. If available, primary calibrators, i.e. plan- 
ets, preferably compact Uranus, were observed at least once 
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Figure 1. The 850 /^m beam map. The contours represent the 
levels -0.3 to 0.9 in steps of 0.2. The data for A2218 were ob- 
tained with different chopping configurations. The beam map was 
constructed by weighting the different chop configurations with 
observing time. Most of the data were obtained with a chop throw 
of 35". The configurations with chop throws of 45" and 90" were 
used for a short time in the observations, which is reflected by 
the few per cent in the beam map. The negative dips at 35" have 
a little tail away from the central peak, which is due to the 45" 
chop throw. 



during an observing shift. The uncertainty in the flux cali- 
bration is approximately 10 % at 850 /im and approximately 
30 % at 450 Atm. 

The on-sky exposure time was 39.7 hours at 850 fim and 
33.4 hours at 450 /im without overhead, i.e. without time 
needed for the chopping, jiggle, etc. The total area surveyed 
is 11.8arcmin'^ at 850 /im and Q.larcmin'^ at 450 /im. 



2.2 Reduction 

The data were reduced using the standard SURF package 
(Jenness & Lightfoot 1998). The reduction procedure is 
described in detail in the catalogue paper for the Leiden- 
SCUBA Lens Survey yKnudson 2004). The maps have been 
regridded with a pixel scale of 1". The full width at half max- 
imum (FWHM) of the beam is 14.3" in the final 850 fim and 
7.5" in the 450 pm. 

Only the beam scale and above is real information re- 
flected. Thus to suppress artificial pixel noise, the maps 
are smoothed to reduce the spatial high frequency noise. 
The 850 /xm and the 450 pm maps were smoothed using a 
2-D Gaussian function with a FWHM of 5" , yielding the 
beams of 15.1" and 9" respectively. Additionally, to match 
the beamsize of the 850 ^m, the 450 /tm map is smoothed 
with a 2-D Gaussian with FWHM of 10" , which results in 
a beam of 12.5" , comparable to the beam size at 850 /im. 
The resulting maps are shown in Figure |21 The sensitiv- 
ity in the final 850 /im map is 0.65mJy/beam in the deep- 
est part of the map and 1.1 mjy /beam averaged across the 
map. This is of equivalent depth to the 850-maps published 
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120021) . The sensitivity in the final 450 /xm 
map is 4.3mJy/beam in the deepest part and 14mJy/beam 
averaged across the map, as determined from the noise sim- 
ulations described below. 

In this paper we concentrate on the 850 /tm map. In gen- 
eral the analysis of 450 /im maps is more difficult because the 
450 fim beam pattern is very sensitive to temperature varia- 
tions on the primary mirror of the JCMT. The beam pattern 
cannot easily be described by a 2D Gaussian function. Fur- 
thermore, both the atmospheric opacity and the uncertainty 
in the fiux calibration are large at 450 nm. This makes ob- 



servations of faint objects difiicult. However, given the good 
quality of the 450 fim data obtained for A2218 we do include 
the 450 fim flux for the 850 /im sources (Table 



2.3 Construction of noise maps 

When a data set is combined and regridded into a map, 
some areas are sampled less than the rest of the map. In the 
A2218 images, these include the areas along the edges, the 
areas where the bolometers are either more noisy or where 
data has been excluded, and areas with less integration time 
outside the overlap between different pointings. It is thus 
imperative to assess the noise as function of position in the 
map to be able to make a reliable analysis of the observed 
flux map. 

When data files for the SCUBA jiggle maps have been 
fully reduced, i.e. all correlated noise has been removed, the 
noise in each individual bolometer can be assumed to be 
independent of other bolometers. The sources we are ob- 
serving are very faint and thus large integration times are 
necessary in order to reach the needed noise levels for de- 
tections. However, the observations are carried out on short 
time scales on the order of an hour, which are the standard 
lengths of our observations, thus the data stream from the 
bolometers are dominated by noise rather than signal. 

In order to determine the position-dependent noise a 
Monte Carlo simulation was used to construct 'empty' (i.e. 
no source signal present) maps representing the noise com- 
ponent in the real data files. For each original data file, the 
statistical properties of each bolometer's data set were cal- 
culated. Then a simulated data set for each bolometer was 
created using a Gaussian random number generator with 
the same properties. This produced an empty map for each 
original data file, and the entire set of these was reduced and 
concatenated using SURF in the same manner as the original 
data files. The resultant single map is thus one realization of 
the empty map simulation procedure. This process was then 
repeated about 500 times. Noise maps were constructed from 
the simulated empty maps by taking all the empty maps and 
calculating the standard deviation at each pixel. The noise 
maps at both 850 /tm and 450 /tm are shown in Fig.|21 where 
it can be seen how the noise is not uniform across the field 
and how the different effects described above appear. 

This algorithm is similar to looking at the scatter of the 
data that goes into each pixel. It does not take into account 
residual 1/f noise and other systematic errors. Comparing 
the two methods suggests the difference is only at the 10% 
level (C. Borys, private communication). 



2.4 Deconvolution of beam map 

As described in section [2.11 the chopping of the secondary 
mirror on either side of the central pointing creates a beam 
pattern with a positive peak and two negative sidelobes. 
This is undesirable because negative sidelobes can hide 
sources or bias their fluxes. Also the negative sidelobes con- 
tain useful signal that can be used in source detection. 
Therefore it is useful to deconvolve the beam pattern out. 
Since we are dealing mostly with isolated p oint sources the 
classical CLEAN algorithm llHogbomlll974f) is ideally suited 
for this purpose. In summary, the classical CLEAN algo- 
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Figure 2. In this mosaic the SCUBA flux maps (top) of A2218 and the corresponding noise maps are shown (bottom). The left side 
represent the 850 fim maps, and the right side the 450 fim maps. Top-left: The 850 fim flux map overlaid with 850 ^m signal-to-noise 
contours of S/N =3,4,5,6,8,12,18,24. The map has been cleaned and sources have been restored using the central positive peak of the 
beam (see Sect. 12.41 . The map has been smoothed with a two-dimensional Gaussian with width 5". The white line indicate the area being 
used in the further analysis. Top-right: The 450 fim flux map overlaid with the A50 fim signal-to-noise contours of S/N =3,4,5,7,10. The 
map has been smoothed with a two-dimensional Gaussian with width 10" . Bottom-left: The 850 fim noise map. This map represents the 
position dependent noise after smoothing with 5". The contours have been overlaid to enhance the contrast and are 0.8 to 3.6mjy/beam 
in steps of 0.4mjy/beam. Darker shades of grey indicate a lower noise. The dashed line indicates the area of the map. The edges are 
clearly less sensitive and the whole edge region of the map is not included in the analysis. The white line indicate the area being used 
in the further analysis. Bottom-right: The 450 fim noise map. This is the noise map for the data smoothed with 10". The contours have 
been overlaid to enhance the contrast and represent 4 to 16mjy/beam in steps of 2mjy/beam. The dashed line indicates the area of the 
map. Comparing left to right, the slightly smaller field of view at 450 /xm is evident. 



rithm is an iterative algorithm which searches for the bright- 
est point, which matches the beam pattern, in the final map. 
The beam pattern is scaled to some fraction of the peak flux 
and subtracted from the map at that position. The whole 
process is then repeated in principle until the residual map 
has an r.m.s. comparable with the noise, though in prac- 
tice until a given flux limit. All the information about sub- 
tracted fluxes and positions (the resultant "delta-functions" ) 
can then be used for source detection and restoration with a 
different beam pattern e.g. the central positive peak of the 
beam pattern. 

In our cleaning algorithm, the (cleaned) data map was 
convolved during each iteration with the beam map and di- 
vided by the associated convolved noise map. The former 
step enhances any present signal in the map whilst the lat- 
ter ensures that no noise peaks (e.g. along the edges) would 
be mistaken for signal. This image is then used to select 
the peaks, but the cleaning itself takes place on the uncon- 
volved data map. We found that subtracting 10 % of the 
flux at each iteration produced acceptable convergence. The 
cleaning was continued until the residual map contained no 
pixels brighter than three times the r.m.s. noise level. Fi- 



nally, the detected sources were convolved with the positive 
peak of the beam map and coadded to the uncleaned resid- 
uals to give a final cleaned map. Using the central positive 
peak to restore the sources in the map preserves the infor- 
mation about the beam shape, which will be important for 
the source extraction. The cleaned 850 fim map, where the 
sources have been restored using the central positive peak of 
the beam, is shown in Figure |21 The 450 fim map has been 
cleaned in the same manner and is also shown in|5| 



3 MEXICAN HAT WAVELET SOURCE 
DETECTION 

Independent of the survey strategy, the source extraction 
algorithm applied to the SCUBA maps must be a robust, 
well-characterized method. Indeed, the noise in SCUBA jig- 
gle maps has a temporal variation and the sources we wish 
to detect typically have at best moderate signal-to-noise ra- 
tio (S/N). The far-infrared emitting regions in high redshift 
sources are expected to be relatively compact and so are 
likely to be unresolved in SCUBA's 15" beam. While Mexi- 
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can Hat Wavelet source detection is the scope of this paper, 
we first summarise a number of techniques have been used 
for locating point sources in previous SCUBA jiggle-maps. 

In the cluster lens surveys, extraction techniques have 
tended to be simple r, as these surveys c over s maller fields 
with fewer sources. |SmaiLg1i_alJ |^22al) and ISmail et alJ 
J2QQ§) used SExtractor l|Bertin fc ArnoutJl99^ . which was 
primarily developed for use in optical images which have 
very different noise ch aracteristics to SCUB A maps. In the 
Hawaiian deep fields, iBareer et al.l (|l99flh simply used a 
signal-to-noise criterion for t heir source selectio n, a tech- 
nique also used for the HDF ijHughes et alj |l933 ). 

For the ongoing SCUBA Half-Degree Extragalactic Sur- 
vey (SHADES) a matched-filter source-extrac tion technique 
has b een demonstrated for a sub-set of data iMortier et alJ 
I200I). 

In their jiggle-map s of the Canada U K D eep SCUBA 
Survey (CUDSS) fields. lEales et al.1 ll200d) and I Webb etlD 
i2003h selected sources from a convolved signal-to-noise 
map, which was constructed by convolving their reduced 
map with a beam map derived from a calibrator, and divid- 
ing by a simulated noise map (see section . An iterative 
deconvolution (CLEANing) routine was then used for deter- 
mination of the position a nd fiux of the dete cted sources. 

In the 8-mJy survey JScott et al.ll2002h . sources were 
extracted by a simultaneous maximum-likelihood fit to the 
fiux densities of all probable peak locations, selected purely 
on the grounds of fiux. A template beam map from a cal- 
ibration observation was centred on every pixel with fiux 
> 3 mjy in a map convolved with a Gaussian filter. The 
height of each potential source was then increased inde- 
pendently until a minimum was found in the statistic 
between this built-up map and the real map. A similar ap- 
proach was used in the extended Hub ble Deep Field (HDF ) 
and flanking fields source extraction llSerieant et alji2003ll . 
In short, in the two methods applied to the CUDSS and to 
the 8-mJy survey data sets the images are convolved with 
the beam. The methods differ in how the convolved images 
are converted to a catalogue. 

3.1 MHW Detection parameters 

The standard mathematics and detection parameters in- 
volved in the Mexican Hat Wavel et (MHW) transform 
were given in iBarnard et al.l (l2004|j . Further details can 
also be found in the original papers describing this work 
JCavon et al.ll200d: IVielva et al.ll2001allbl) . 

Point source detection is primarily controlled by setting 
numerical requirements for two parameters, namely an opti- 
mal scale, Ropt, and an wavelet coefficient value, wf[Ropt)- 
We performed initial, simple simulations, using both real 
data and the simulated noise maps from the previous sec- 
tion, to understand the response of the MHW routines to 
SCUBA jiggle-maps, since these are quite different in sev- 
eral respects to the simulated Planck data sets on which the 
technique was previously tested. The procedure involved in 
selecting sources for the final detection list is as follows: 

(i) Firstly the optimal scale, Ropt, is calculated by the 
MHW software for each input map. This involves iterating 
through small changes in the value of R around the point 
source scale, 9, until the maximum amplification is found 



for the map. The final value will be depende nt upon the 
measured impact of the noise on varying scales dVielva et alJ 
l2001b^ — noise with a characteristic scale a little larger than 
9, for instance, can be most strongly counteracted by an 
optimal scale -Ropt slightly smaller than 9. 

(ii) Point source candidates are then selected at posi- 
tions with wavelet coefficient values Wf{Ropt) ^ 2(Tii,„ (Ropt), 
where ajn„ (Ropt), where a^^ (R) is the dispersion of the noise 
field in wavelet space. The value of 2 was suggested by our 
early simulations as a value which allowed all real sources 
to pass through to an initial catalogue. 

(iii) For each candidate, the 'experimental' Wf{R) is com- 
pared to the theoretical variation expected with R, as a fur- 
ther check on the source's shape. A value of is calculated 
between the expected and experimental results and the sec- 
ond parametrised requirement is therefore that ^ 4, i.e. 
that the region surrounding the identified peak has the char- 
acteristics of a Gaussian point source. 

Since, as discussed in Section f2. 31 the noise level in the 
map varies widely, the input map to the MHW software 
used was a signal-to-noise map, created by dividing the final, 
smoothed data map by the final, smoothed noise map. This 
normalises the input map with respect to the noise. This is 
especially useful for this technique as the properties of the 
whole map are used to determine detection parameters, and 
so high noise values can distort these values unfairly. 

Despite this normalisation, detections close to the edge 
of the map, where the noise levels are at their highest, were 
considered to still be unreliable and such detections in a 1.5- 
beam region around the edge are not retained in the final 
catal ogue. As discussed in previous papers fe.g.. llvison et alJ 
|2_002f) the number of spurious detections in SCUBA maps is 
larger along the edges than elsewhere in the maps. 

To assess the chances of obtaining false positive detec- 
tions in the S/N map, we performed the MHW source ex- 
traction upon empty S/N maps, created by dividing the 
empty Monte Carlo maps by the final average noise map. In 
100 such realisations we found only two sources with S/N > 
3. We thus conclude that, when restricting the sources to 
those which both meet the standard detection parameter re- 
quirements outlined and which have S/N > 3 in real space, 
the probabililty of spurious detections is very low. A detailed 
analysis of this is performed for the full Leiden-SCUBA Lens 
Survey (Knudsen, 2004; Knudsen et al., in prep.). 

3.2 Comparison of MHW and CLEAN source 
detection 

We compare the MHW source detection done on the un- 
cleaned map and the cleaned map. x,y pixel positions rela- 
tive to the centre of the map and fiuxes for detections are 
fisted in Table[T] The first nine sources are numbered accord- 
ing to their number in the final catalogue as given in Table|21 
The positions agree very well for all nine sources, though for 
source 7 there has been a four pixel shift. Source 7 is close 
to the negative sidelopes of the bright central sources. In 
the cleaned map its wavelets parameters are much improved 
and its profile provides a better match to the theoretical 
expectation as estimated by the MHW algorithm. Source 7 
is the reason why we CLEANed the SCUBA map for this 
project, i.e., to confirm whether the source was an artifact 
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Table 1. MHW and CLEAN detections 



MHW CLEAN 
uncleaned cleaned 



src 


Ax 


A?/ 


/ 


Ax 




/ 


Ax 


Ay 


/ 


1 


78 


-40 


10.8 


79 


-40 


10.4 


78 


-41 


9.0 


2 


20 


-17 


9.0 


20 


-17 


8.7 


18 


-17 


9.2 


3 


-1 


2 


17.1 


-0 


1 


16.1 


-1 





22.3 


4 


-6 


15 


11.7 


-6 


14 


12.8 


-8 


14 


11.3 


5 


-6 


-34 


5.0 


-6 


-34 


3.1 


-7 


-34 


3.2 


6 


-8 


36 


11.0 


-8 


36 


11.3 


-11 


36 


13.9 


7 


-54 


30 


4.3 


-51 


31 


3.3 








8 


-69 


35 


5.4 


-69 


35 


5.2 


-68 


32 


3.8 


9 


-75 


10 


4.6 


-74 


10 


4.8 


-75 


10 


4.4 


10 


45 


13 


4.4 














11 


59 


-17 


4.0 














12 


113 


-1 


9.3 














13 














-30 


68 


2.0 


14 














-10 


69 


2.3 


15 














2 


29 


2.4 


16 














30 


50 


2.6 



Note: Ax and Ay give the distance in pixels from the map cen- 
tre at a,(5 = 16'^35™54.22=,-|-66°12'24", and the pixel scale is 1 
arcsec/pixel. The flux, /, is in units mjy. 

due to the negative sidelopes of the beam. For the other 
sources the wavelets parameters were improved only by a 
few percent when using the cleaned map. In the uncleaned 
map three additional sources were detected. Source 12 falls 
within the edge region, which has been trimmed off. Sources 
10 and 11 do not produce a good match to the theoretical 
wavelets expectation. As these sources were not re-detected 
in the cleaned map, we consider them spurious. 

Furthermore, we compare the MHW source detections 
with the detections from the CLEAN algorithm (also shown 
in Table Ql. The CLEAN algorithm detects twelve sources, 
of which eight are in good agreement with those detected by 
the MHW algorithm. CLEAN yields larger flux estimates for 
sources 3 and 6. The positions deviate by about 1-2 pixels 
from the MHW position, but in general there is good agree- 
ment between the two methods. Source 7 is not detected 
by CLEAN. The four additional sources (13-16) are faint, 
~ 2 mJy, and remain undetected by MHW. The reason for 
this could be found in how the two methods depend on the 
zero-point of the map: If large scale variations across the 
map elevate the zero-point of the map, this will also have an 
impact on the S/N value so that these appear systematically 
off-set by a small value. For the CLEAN routine, which se- 
lects the sources by S/N > 3, also positions with artificially 
increased S/N ratio will be selected. The MHW algorithm, 
which searches for sources of a given scale placed on top 
of a background, such large scale variations will have no or 
only very small infiuence on the detection of sources, and 
thus MHW remains essentially insensitive to the zero-point 
of the map. For future large surveys, such as those planned 
to be carried out with SCUBA-2, large scale variations can 
pose a problem. MHW may offer a good solution for source 
extraction from such maps. 

Finally, in order to check the MHW source detection, 
in order to check for possible sources hidden in the wings 
of bright sources, we subtracted the detected sources from 
the map, using a template for the central beam scaled to 



the fluxes of the individual sources. This residual map is 
not the residual map from the CLEAN routine, as it still 
contains the sources which were detected by CLEAN but 
not by MHW. Then, we performed source extraction with 
MHW on the residual image. If there is a source closer than a 
beam to the other sources, which has remained undetected, 
then we expect to detect it in the residual map. Four such 
sources are detected. Three of these are directly associated 
with bright sources and are likely artifacts resulting from an 
imperfect match of the bright source and the scaled beam 
subtracted from the map. The fourth source may be real, but 
is outside of the 7.7arcmin^ trimmed region near the SW 
edge and thus no included in the analysis. So in conclusion 
we believe that our MHW catalogue is robust and does not 
miss any sources to the detection limit of the data. 

3.3 Testing MHW with simulations 

Using the empty Monte Carlo maps we investigated the de- 
tection accuracy of the MHW algorit hm in a similar fashio n 
as was performed for the scan-maps in lBarnard et alJi2004l) . 
One point source was added at a time to the empty maps 
using the central peak of the beam map shown in Figure 
as a template source. This was repeated 400 times for each 
flux step. The position was randomly chosen, with a uni- 
form distribution across ten maps, so that the whole map 
area was well sampled. The input flux of the point source 
was then increased and the standard MHW detection algo- 
rithm was run on each map. Additionally, we also perform 
two simulations for the CLEAN routine concerning sources 
close to the edge and overlapping sources. 

3.3.1 Position accuracy 

In Figure |21 we plot the standard deviation of the differ- 
ence between input and detection position. The standard 
deviation decreases with increasing S/N and has a value 
about 2.5" at S/N ~3 decreasing towards 0.5" for large 
S/N . We adopt these as the uncertainty on the determined 
position introduced by the extraction algorithm. This added 
in quadrature with the pointing uncertainty of the JCMT 
(~2") gives the error on the determined position. 

3.3.2 Flux accuracy 

In Figure |21 we plot the relative average flux difference, 
{Sout — Sin)/ Sin, and the standard deviation as function 
of S/N and of flux. For S/N > 7 and also flux > 7 mJy the 
detected flux on average is overestimated by one to two per- 
cent. Eddington bias is seen through the average overestima- 
tion of the sources with fainter fluxes and lower S/N . This 
is also seen in oth e r SCUBA surveys f e.g..|Scott et a l]|2002|: 
IWebb et all l2003l : iBorvs et al] l2003l: iBarnard et alJl2004^ . 
We note that the calibration uncertainty is larger than the 
Eddington bias. The standard deviation on the average rel- 
ative flux difference decreases with S/N, from about 25% at 
S/N = 3 to 5-7% at S/N > 15. As uncertainty on the de- 
termined flux we use the standard deviation of the relative 
average flux different, which is added in quadrature with 
the calibration error from the SCUBA reduction to give the 
total error on flux density of the individual sources. 
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Figure 3. Upper panels: The relative difference between the input and detected flux as function of S/N (left) and flux (right) based on 
the simulations. The flUed dots represent the average in each S/N or flux bin, while the error bars represent the standard deviation in 
the same bin. Lower left panel: The standard deviation (•) of the difference between input and detection position as a function of S/N 
based on the simulations. Lower right panel: The completeness as function of input flux based on the simulations. The 850 fim map is 
80% complete at 3.8 mjy and 50% complete at 2.7 mjy. For comparison, 3mjy corresponds to S/N ~ 3 in about two-thirds of the map. 



3.3.3 Completeness 

The completeness as function of input flux was determined 
from the same set of simulations. The fraction of sources 
detected at each flux level is calculated. The result is plotted 
in Figure 13 The observations are 80% complete at 3.8 mJy 
and 50% complete at 2.7mJy. 

3.3.4 Sources near the edge 

For sources closer to the edge one of the two negative side- 
lopes will be outside of the map. This can complicate the 
source detection in CLEAN, where a matching to the full 
beam pattern is done. In simulating 400 point sources, one 
at a time, each of 20 mJy, we found that 68 sources, i.e. 
17%, were not detected, while similar sources away from the 
edge were all detected. Additionally, the position and flux 
had larger uncertainties compared to those away from the 
edge. 

3.3.5 Overlapping sources 

Given the large width of the 850 /im beam, blending of 
sources is not an uncommon occurance. For example in the 
A2218 two of the nine sources are separated by approxi- 
mately one beam width. We simulated source blending by 
adding two sources with a separation between and 30 arc- 
sec onto the MC maps. The flux ratio of the two sources 
was taken to be 1:1, 1:2 and 1:4, although we flnd that 
different flux ratios do not significantly alter the result in 



either case. MHW detects individual sources separated by 
more than one beam width, while CLEAN is able to detect 
sources with separations > 12". However, the flux and posi- 
tions accuracies using MHW are on average approximately 
5% better than those of CLEAN. 



4 ANALYSIS OF ABELL 2218 SUB-MM MAPS 

4.1 The sub-mm source catalogue 

We now present the detailed analysis of the Abell 2218 
SCUBA maps using the MHW technique. We focus on a 
total area of useful data of 7.7arcmin^ (see flgure|5J. Nine 
sources were detected by MHW using the detection param- 
eters described in Section 3.1. Six of these sources have de- 
tectable 450 /im emission. The detected sources are listed 
in Table |21 The sources are named using the preflx SMM 
followed by their J2000 coordinates. For the 850 fim param- 
eters we give the S/N in real space, the combined flux and 
flux uncertainty and position uncertainty. For the 450 fim, 
where no emission is detected, we give 3cr upper limits. The 
uncertainty on the 450 fim flux is 30%, which is derived from 
the calibration of the absolute flux. In Figure 2] we show 
the 850 /im S/N map, and the S/N contours overlaid on an 
optical image for orientation. The differences seen in flux, 
in particular at 450 /im, between these results and those of 
iKneib et all ll2004a'l can be assigned to uncertainties in the 
calibration. More data has been in cluded in th e combi ned 
data set presented here than in the lKneib et al.l (1200481) . 
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Table 2. Sources detected in the SCUBA maps of A2218. / is flux; S/N is the signal-to-noise of the source in real space as detected 
using MHW; cr gives the uncertainties on the flux and position at 850 and 450 /^m respectively. The flux uncertainty is that found for 
the appropriate S/N bin added in quadrature with the calibration uncertainty, as described in section r3.3.2l The positional uncertainty 
does not include confusion errors. 



NAME 




pos 


/SSO ± O-850 


SIN 


/45O ± O"450 


SIN 






/85o/M 






U 


mJy 




mJy 


mJy 






mJy 


SMM J163541. 2+661144 


(1) 


2.9 


10.4±1.4 


7.5 


< 59.4 




3.188 


1.7 


6.0 


SMM J163550.9+661207* 


(2) 


2.7 


8.7±1.1 


11.5 


39.2±11.8 


7.2 


2.516 


45* 


0.8 


SMM J163554.2+661225* 


(3) 


2.6 


16.1±1.6 


21.7 


89.7±26.9 


19.0 


2.516 


45* 


0.8 


SMM J163555.2+661238* 


(4) 


2.6 


12.8±1.5 


16.9 


60.1±18.3 


12.1 


2.516 


45* 


0.8 


SMM J163555.2+661150 


(5) 


3.5 


3.1±0.7 


3.8 


29.1±8.7 


6.1 


1.034 


7.1 


0.4 


SMM J163555.5+661300 


(6) 


2.6 


11.3±1.3 


15.8 


18.0±5.4 


3.0 


4.048 


4.2 


2.7 


SMM J163602.6+661255 


(7) 


3.5 


2.8±0.6 


3.5 


< 21.6 






1.8 


1.6 


SMM J163605.6+661259 


(8) 


3.3 


5.2±0.9 


4.9 


< 25.5 






1.8 


3.3 


SMM J163606.5+661234 


(9) 


3.3 


4.8±0.8 


4.6 


30.1±9.3 


3.1 




1.6 


2.9 



* SMM J163550.9+661207, SMM J16 3554.2+661225 and SMM J163555. 2-1-661238 have been identified as a multiply-imaged galaxy, 
which is denoted SMM J16359+6612 iKneib et al.ll2004a^ . 



4.2 Gravitational Lensing Magnification 



We e xploit the detailed mass model of A2218 ijKneib et alJ 
Il996|) updated to include the triple sub-mm image 
jKneib et aDl2004ari and the high redshift multiple images 
at 2 = 5.56 (Ellis et al., 2001) and 2 ~ 7 (Kneib et al 2004b). 
In Figure |S| we show the area as function of magnification 
for source planes at the redshifts z = \ and 4. In 1.3 arcmin^ 
the flux magnification factors are > 2. Furthermore, in Fig- 
ure|K|we show the area as function of sensitivity both in the 
image plane and in the source plane, again for source planes 
at different redshifts, 2 = 1 and 4. For high redshift, > 1, the 
lensing magnification is only weakly depending on redshift. 
For the sources without known redshift, assuming a redshift 
> 1 will not have a strong impact on the derived properties. 
For z > 2, the effective area surveyed (source plane) is about 
2.7arcmin^ which corresponds to an average magnification 
factor of 2.8. 

The redshift is known for six sources in the field: three 
of these are a multiply-imaged object at redshift z = 2.516 
(Kneib et al., 2004a), while the other three are sin gle-imaged 
sources at redshift z = 1.034 jPello et alJll993) . and z = 
3.188 and 4.048 (Paper II). For the sources with unknown 
redshift we assume z — 2.5, based on the median redshift 
for br i ght S CUBA sources determined by IChapman et alJ 
l|2003l l2005tl . The magnification factors of the individual 
singly lensed sources range between 1.6 and 7.1. For the 
multiply-imaged source at z = 2.516 the magnification fac- 
tors for individual images are 9, 14 and 22, i.e. a total mag- 
nification of 45 for the three images combined (Kneib et al., 
2004). In Table|5|we list the magnification factors, /i, and the 
lensing corrected fiuxes, /sso/m for the individual galaxies. 

Note that three (five, including all images of 
SMMJ16359-I-6612) of the nine detected sources have un- 
lensed fiuxes fainter than or close to the blank field confu- 
sion limit. These sources would most likely not be detected 
in blank field surveys and thus not have been accessible to us 
without the use of gravitational lensing or a larger telescope 
with sensitive instruments. 



4.3 Confusion 

Confusion noise is caused by unresolved faint sources in the 
field. In the background of deep SCUBA maps the instru- 
mental noise and the confusion noise from a fainter submm 
population are of approximately equal magnitude. Under- 
standing the confusion in a map is essential to be able to es- 
timate the relative importance of uncertainties due to blend- 
ing of sources. 

For the cluster fields the confusion limit is affected by 
the gravitational lensing. The gravitational lensing magnifies 
the region seen behind the cluster, hence the source plane 
is smaller than the image plane. The number of beams is 
conserved between the image plane and the source plane, i.e. 
the size of the beam scales with the magnification. This is 
why it is at all possible to observe the fainter sources, which 
have a higher surface density than the brighter sources. 

The number counts in the lensed case can be written 
as Niens = {Nolfi){Slfj,)"°' = NbiankfJ'"''^, whcre is the 
gravitational lensing magnification, assuming a blank-field 
number count power law where A'^(> 5") = NqS'" and S 
is the observed fiux. The average magnification for a field 
can be found as the ratio of the area in the image plane 
and the area in the source plane. We use the rule of thumb, 
that the confusion limit in ima ging is defined as one source 
per 30 beams (e.g. iHogdliool . The confusion limit in the 
lensed case can thus be written as Sconf ~ {SOQNo^j,^^")^^" 
(illustrated in Figure |KJ , where Q is the solid angle of the 
beam. As the lensing magnification varies across the field, 
this means that the actual confusion limit is deeper in highly 
magnified regions close to the caustics than the estimated 
average lensed confusion limit. 

In case of the A2218 850 fim map, the confusion limit in 
the source plane based on this simple estimate is ~ 1.1 mJy 
on average across the whole field. In the simplified esti- 
mate of the confusion limit presented here we assumed that 
the number counts are described by a single power law. 
There are good indications that the number counts are de- 
scribed by a double power law or another function with a 
(gradual) turn-over llScott fc Whitj|l999l: IScott et aljbool 
iBorvs et alJl2003l: iKnudse J2OO4I: Knudsen et al., in prep.). 
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Figure 4. Upper panel: A2218 850 fim S/N map where the 
nine detected sources each are indicated with a circle with diam- 
eter of 15". The number assigned to each source corresponds to 
their order in Table 1^ The white line indicates the area used for 
the analysis. Lower panel: The 850 /xm S/N = 3, 4, 5, 7, 10, 15 
contou rs overlaid on the HST ACS z-band image iKneib et alJ 
l2004bft . 



In the Leiden-SCUBA Lens Survey ((Knudse^lool) , the 850 
number counts are probed over almost two decades of flux 
(0.1-20 mjy) and improving the statistics in particular on 
the faint end. Including double power law number counts in 
such a calculation will work in a favourable direction and 
the confusion limit in the source plane will be lower. 

4.4 850 fim source counts in A2218 

In Fig. 1^1 we plot the number counts for sources detected in 
A2218. The number counts are calculated for the flux levels 
corresponding to the flux of each of the sources. Because 
of the large variation in the sensitivity across the map, it 
is not possible to directly determine the number counts as 
cumulative number counts. The complic ations arising from 
this has been disc ussed in previously (e.g.. lBorvs et al.l2003l : 
IWebb et alJl2003ll . We here apply a simple method for esti- 
mating the number counts: At each flux level, S, only the 
area where S > 3ct are considered. The number of sources 
within that area are divived by the area. We use the sen- 
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Figures. Upper panel: The area of the source plane as function 
of magnification for redshifts z = I (dashed) and z = 4 (dash- 
dot). The difference between the different redshifts is relatively 
small. Middle panel: The area of the source plane as function of 
3(T sensitivity for redshifts z = 1 (dashed) and z = 4 (dash-dot). 
The source plane a.t z = 2.5 is about 2.8 times smaller than the 
observed area of the image plane. Lower panel: The confusion 
flux limit as function of magnification, as described in subsection 



sitivity map for the source plane at redshift z — 2.5 and 
note that using that of sources planes at other redshifts does 
not make a noticable difference. By using the source plane 
sensitivity map, we also take into account the effect that 
lensing has on the area. The error bars are Poisson statistics 
jGehrelslll986D . The number counts are in reasonable agree- 
ment with the results from other surveys. Differences be- 
tween surveys can be assigned to cosmic variance and small 
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number statistics. The 850 number counts at the level 
of ~ ImJy is ^ 2arcmin~^. Additionally, we determine an 
upper limit for the number counts at Sgso = 0.2 mjy and 
0.1 mjy of 44arcmin~^ and 175 arcmin"^ respectively. The 
450 At m number counts p r obe deeper than p reviously pub- 
lished ISmajTeFaO ll2002l) : lBorvs et al.l ^|2003^ . They appear 
to be slightly higher, which can be assigned to sample vari- 
ance, but agree within the error bars. 

The fraction of the submm extragalactic background 
light (EBL) resolved in the SCUBA map can be constrained 
through the sum of the fluxes. The total flux relative to the 
area is (3 5±1.6) Jydeg-^ Thi s is ~ 80% of the EBL as deter- 
mined bv lFixsen et al. (19 9^, and co mparable to the EBL 
value determine bv iPuget et al.l lll99d) . This al most double 
the to tal flux value of 17.3 Jydeg"^ found bv lCowie et alJ 
l|2002l) . Due to cosmic variance and the clustering of SCUBA 
sources, this value is expected to vary between fields. How- 
ever, the fact that we resolve such a significant fraction of 
the EBL, is a good demonstration of the power of gravita- 
tional lensing in combination with the SCUBA observations, 
as it allows us to probe the faint submm galaxy population, 
which is the dominant contributor to the submm EBL. 



5 CONCLUSIONS AND OUTLOOOK 

We have performed a thorough analysis of the Mexican Hat 
Wavelet algorithm (MHW) as a source extraction tool ap- 
phed to SCUBA jiggle-maps. The analysis was done us- 
ing the deep SCUBA maps of the galaxy cluster A2218, 
through source extraction with MHW on the real data and 
through extensive simulations. We found that MHW is a 
stable method for source extraction at low signal-to-noise 
{S/N > 3). We conclude that MHW is an algorithm suitable 
for source extraction from SCUBA jiggle maps and is a pow- 
erful tool for studying the faint submillimetre sources and 
thereby the faint end of the submillimetre number counts. 
MHW has potential as source extraction algorithm for data 
taken with the future SCUBA-2 instrument. SCUBA-2 will 
make "total power" maps, which should be free of chopping 
artifacts, hence MHW will have an immediate advantage 
over the other techniques currently in use. 

The SCUBA map we have obtained for A2218 is one 
of the deepest submillimetre maps ever taken. Covering a 
large area of the cluster, we have been able to survey the re- 
gion where strongly lensed background sources are present. 
In the analysis of A2218, nine sources were detected in the 
850 /^m SCUBA map, the largest number detected in a sin- 
gle ultra-deep SCUBA map. Six of these sources were also 
detected at 450 /^m. Correcting for the gravitational lensing 
by the galaxy cluster, three sources have intrinsic 850 /im 
fluxes below the blank field confusion limit, and three with 
fluxes comparable to the blank field confusion limit. It is 
the presence of strong gravitational lensing, which pushes 
the confusion limit to fainter fluxes levels, which has made 
it possible to detect the faint submillimetre source with 
< 2 mJy. For this field we determine the 850 /xm number 
counts to Sgso = 0.4, and place upper limits at Sgso = 0.1 
and 0.2 mJy. Additionally, we also estimate the 450 /im num- 
ber counts down to 5450 = 4mJy. The identification of the 
individual sources will follow in Paper II. 

ALMA will in the future allow for yet deeper observa- 



tions and thus a continuation of the number counts to even 
fainter fluxes. If ALMA will detect a number density of ~ 50 
galaxies/arcmin^, this means we are reaching the number 
density of galaxies detected in the optical/near-infrared and 
thus will allow for a better understanding of the connec- 
tion between dust and stars in the Universe as function of 
time/redshift. 
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